Three-dimensional Casimir force between absorbing multilayer dielectrics* 
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Recently the influence of dielectric and geometrical properties on the Casimir force between dis- 
persing and absorbing multilayered plates in the zero-temperature limit has been studied within a 
ID quantization scheme for the electromagnetic field in the presence of causal media [R. Esquivel- 
Sirvent, C. Villarreal, and G.H. Cocoletzi, Phys. Rev. Lett. 64, 052108 (2001)]. In the present 
paper a rigorous 3D analysis is given, which shows that for complex heterostructures the ID theory 
only roughly refiects the dependence of the Casimir force on the plate separation in general. Further, 
an extension of the very recently derived formula for the Casimir force at zero temperature [M.S. 
Tomas, Phys. Rev. A 66, 052103 (2002)] to finite temperatures is given, and analytical expressions 
for specific distance laws in the zero-temperature limit are derived. In particular, it is shown that 
the Casimir force between two single-slab plates behaves asymptotically like instead of d~* (d, 
plate separation). 

PACS numbers: 03.70.-|-k, 12.20.Ds, 42.60.Da 



I. INTRODUCTION 

The Casimir effect has drawn attention to itself over 
decades (for a review and a collection of references, see 
Q Q)' Various concepts and calculational techniques 
have been developed, but only few of them have turned 
out to be capable of dealing with realistic (at least 
Kramers-Kronig consistent) dielectric bodies. The famil- 
iar 'mode summation' method of deriving Casimir forces 
as employed by Casimir himself 0], for instance, suffers 
from the obvious fact that there are no modes to ex- 
pand the (macroscopic) field operators when absorbing 
bodies are involved rather than idealized boundaries. A 
generalized mode decomposition (with the mode func- 
tions being orthogonal with respect to an appropriate 
'weighted' inner product) is appropriate only for (narrow) 
bandwidth-limited fields when dispersion and absorption 
may be neglected j4j. By contrast, all frequencies should 
be included in the calculation of the Casimir force. More- 
over, the annoying infinities one has to cope with when 
using modes always demand the application of regular- 
ization techniques P]. To tackle realistic problems, one 
can proceed in essentially three ways: 

1. The electromagnetic field and the material bodies 
are treated macroscopically, and explicit field quan- 
tization is avoided by invoking statistical thermo- 
dynamics to write down the field correlation func- 
tions that are needed in the Maxwell stress tensor. 

2. The electromagnetic field and the material bodies 
are quantized on a microscopic level. The bodies 
are described by appropriate model systems that 
feature dissipation of energy into a heat bath, and 
the coupled field-matter equations are tried to be 
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solved. In general, such a microscopic approach 
to the problem is hardly feasible without rather 
sweeping assumptions about the microscopic pro- 
cesses involved. 

3. The presence of dielectric bodies is described by 
means of a spatially varying permittivity that is a 
complex function of frequency, without addressing 
a specific microscopic model. For given permittiv- 
ity, the macroscopic, medium-assisted electromag- 
netic field is quantized, by using a source-quantity 
representation of the field in terms of the classical 
Green tensor and an infinite set of appropriately 
chosen bosonic basic fields. 

The first method was introduced by Lifshitz He 
obtained the force acting on two (semi-infinite) absorb- 
ing dielectric walls by calculating the Maxwell stress ten- 
sor in the region between the walls from field correla- 
tion functions, which he could evaluate by employing a 
dissipation-fiuctuation relation for the 'random electric 
field' that acts as a Langevin noise source and balances 
the effect of absorption in the Maxwell equations. Al- 
though Lifshitz' calculation is not really quantum, he was 
fully aware of the fact that the "... zero point vibrations 
of the radiation field" cause the effect at zero tempera- 
ture. Lifshitz' results were later re-derived by Schwinger 
et al. using source theory to circumvent quantization p| . 
In a recent paper, Matloob 01 has essentially followed 
Lifshitz, by postulating field correlation functions, with- 
out explicit field quantization. 

An instructive treatment of the Casimir effect on the 
basis of the second method was given by Kupiszewska |1] , 
using a harmonic-oscillator model for the matter. The 
effect of the heat bath accounting for dissipation was 
properly subsumed within a Markovian damping term 
together with a Langevin noise source corresponding to 
the random field in Lifshitz' approach. Unfortunately, 
the calculations were carried out only for one-dimensional 
(single-slab) systems. Recently, this theory has been used 
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to study the Casimir force between multilayered plates 
i 

Here, we base the calculations on the third method 
recently used by Tomas [lOj . It needs no external in- 
put other than the (phenomenologically given) permit- 
tivities of the bodies. Apart from dropping the bulk 
part of the Green tensor (giving rise to unobservable bulk 
stress only), no regularization is needed. Removal of the 
bulk stress, which is sometimes referred to as "remov- 
ing the Minkowski contribution" in the literature, is in 
fact necessary in all approaches, even in the framework 
of Schwinger's clever source theory |^. 

Let us briefly comment on two further techniques that 
have been employed to calculate Casimir forces. In the 
surface-mode approach originally used by van Kampen 
et al. to obtain Casimir forces in the non-retarded 
limit, the calculations are based on a complete set of so- 
lutions to Laplace's equation. Summing, by means of a 
complex contour integral, the zero-point energy assigned 
to each of these solutions then yields the potential energy 
of the force. Later, the method was extended to include 
also retardation |l2l |. However, in this case all normal 
modes and not only the surface modes should be consid- 
ered. Though the method is intrinsically based on mode 
expansion and hence on real permittivities, correct re- 
sults can also be found for absorbing material if at some 
stage of the calculation the permittivities are allowed to 
become complex. In the scattering approach [13>] , which 
employs elastic (i. e. unitary) scattering theory, the sum 
of the mode frequencies is represented as an integral con- 
taining scattering phase shifts or related quantities (see 
also jl|), and the integration is performed in the com- 
plex plane. Again, if complex permittivities are plugged 
into the resulting expressions (derived for non-absorptive 
materials), correct results can be obtained. Note that 
within the framework of both methods the sum over the 
relevant mode frequencies appears as the singularity con- 
tributions to certain complex contour integrals, which are 
more generally valid than the initial expressions based on 
the concept of normal modes. 

In this paper, we first derive a formula for the Casimir 
force between multilayer dielectric plates at finite tem- 
peratures, which may be regarded as being a gener- 
alization of the Lifshitz formula to multilayer sys- 
tems and an extension to finite temperatures of the zero- 
temperature result derived by Tomas In particular 
for one-dimensional systems, the formula reduces to an 
expression of the type derived by Kupiszewska @ . Trans- 
forming the formula to a form that is very suitable for 
further analytical and numerical evaluation, we study the 
dependence of the Casimir force on various system pa- 
rameters (such as the stacking order and the frequency 
response of the permittivity) in detail. The numerical 
calculations are performed for single-resonance dielectric 
matter of Drude-Lorentz type. It is well known that the 
Casimir force between semi-infinite dielectric walls be- 
comes proportional to d"* for large wall separation d. We 
show that is also the 'generic' behavior for layered walls 



in general. However, yet different types of long-distance 
laws are also possible. In particular, for single-slab walls 
of finite thickness, the asymptotic d^"* law changes to a 
d^^ law. We further show that the case of small wall sep- 
aration d can be treated within Lifshitz' approximations 
i 

The paper is organized as follows. In Section^the for- 
malism is outlined and the basic formula for the Casimir 
force is given. Section UTTl presents analytical results for 
large and small wall separation, and Section llVl is devoted 
to the numerical results. A summary is given in SectionlVl 
followed by five appendices. Appendix 1X1 provides some 
basic relations needed for the finite-temperature calcula- 
tion. Useful recurrence relations for the generalized re- 
flection coefficients are given in Appendix^ The prob- 
lem of switching to imaginary frequencies in the basic 
integral expression for the Casimir force is addressed in 
Appendix [CI In Appendix ^ the formalism is applied 
(for comparison) to one-dimensional systems, and in Ap- 
pendix^ a special integral is evaluated. 



II. CASIMIR FORCE 

A. Quantization scheme 

Let E(r), D(r), B(r), and H(r) be the medium- 
assisted (macroscopic) electromagnetic field operators 
such that 



E(r) = / du;E(r,w) -l-H.c, 



(1) 



and D(r), B(r), and H(r) accordingly. Within the frame- 
work of the quantization scheme given in Ref. |l4j . the 
operators E(r, w), D(r,w), B(r, w), and H(r, w) obey 
Maxwell's equations 



V X E(r, iu) - iLuB{r, lu) = 0, 

VB(r,c.) -0, 

V X H(r, uj) + iwD(r, lu) = 0, 

VD(r,w) = 0, 



(2) 
(3) 
(4) 
(5) 



where (for non-magnetic) linear media the constitutive 
relations read as 

D(r, uj) = eoe(r, c^)E(r, lu) + P^{r, uj), (6) 
H(r,u;)=Mo"iB(r,c.). (7) 

Here, the complex permittivity 

e{r,uj)^e'{r,uj)+ie"{r,Lu) (8) 

satisfies the Kramers-Kronig relations, and the noise po- 
larization 

P^{r, uj) = i^/heoe"{r,uj)/TT{{r, cj) (9) 
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associated with material absorption is expressed in terms 
of bosonic fields f(r,w), 

[/,(r,c.),/,(r',L.')]=0, (10) 
[Mr,^)j]{r',uj')]^S,,6{r-r')5{u;~u;'), (11) 

which play the role of the dynamical variables of the com- 
bined field-matter system consisting of the electromag- 
netic field, the medium polarization, and the heat bath 
accounting for absorption. From Eqs. |2Jl - lEJl it then 
follows that 

B(r, uj) = {iLuy^V X E(r, uj), (12) 
D(r, uj) = inauJ^y^V x V x E(r, w), (13) 

where 



E(r,c^) - 




')f(r',^), 



(14) 

with G(r, r' ,uj) being the (classical) Green tensor, which 
is determined by the equation 

V X V X G(r, r', uj) - ^ e(r, uj) G(r, r', uj) = ^(r, r') 

(15) 

[(5(r,r'), tensorial J-function] together with 'outgoing' 
boundary conditions. In particular, the Green tensor sat- 
isfies the relations 14] 

G,,(r,r',c^) = Gj,(r',r,u;) (16) 

and 

— / d^se"{s,uj)G{r,s,uj)G*{s,r',uj)^^[G{r,r',uj)] 

(17) 

(3, imaginary part). The above given equations do not 
refer to a specific picture. In particular, in the Heisenberg 
picture the operators f (r, uj) simply carry an exponential 
time dependence, 

f(r,w) f(r,t^)e-'"*, (18) 

according to the Hamiltonian of the (overall) system 

H ^ dujhuj J cfrf''{r,uj)f{r,uj). (19) 

For the multilayer structure sketched in Fig. ^ and con- 
sidered throughout the paper, 

£(r, uj) = ei{ui) for zi < z < zi+i = zi + di (20) 

{I = 0, . . . ,n, Zq = — oo, Zn+i = +oo), the Green tensor 
Gi(r,r',Ci;) in the Ith layer can be decomposed as 

G,(r,r',c.) = Gl^""^(r,r',c.) + Gr*(r, r', c^), (21) 

where Gj'""^ is the solution emerging from a point- 
like source current placed in the /th layer without any 




FIG. 1: Scheme of the multilayer structure. 

boundaries, and Gf^'^'^ solves the homogeneous version of 
Eq. H15() so as to make the full Green tensor obey the cor- 
rect boundary conditions at the surfaces of discontinuity. 
Since there is no Casimir effect in homogeneous space, 
the bulk part can safely be dropped (and actually must 
be dropped) in the stress tensor. The remaining scatter- 
ing part contains the geometrical data of the problem and 
is continuous at r = r', whereas the bulk part is rather 
singular there. The scattering part can be constructed 
as shown in Ref. Here, only the case where both 

spatial arguments are in the same layer is of interest. 

B. Stress tensor 

In order to determine the Casimir force from the stress 
tensor, we first calculate 

T(r,r',t) = Ti(r,r',i)+ T2(r,r',i) 

-iJTr[ri(r,r',i)+ T2(r,r',t)], (22) 

where 

ri(r,r',i) = (D(r,f)®E(r',i)) (23) 

and 

r2(r,r',t) = (B(r,t)®H(r',i)) (24) 

(rT^r'). Here, the electromagnetic- field operators are 
thought of as being expressed in terms of the fundamen- 
tal fields f (r, uj) as outlined in Section III Al To specify 
the quantum state, we assume thermal equilibrium. 

1. Basic equation 

For finite temperatures T, we may employ the canon- 
ical density operator 

p = Z-y-"'^'^-^\ (25) 
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where 

Z = Tre-^/('^-^). (26) 
After some calculation we derive (Appendix ^ 
ri(r,r') =Tr[pD(r,t)(g)E(r',t)] 

^-j^ da;coth(^2fc^j!^S[e(r,c.)G(r,r',a;)], (27) 
and 

T2(r,r') = Tr[/5B(r,i)®H(r',t)] 



dcj coth 



V X 3[G'(r,r',t^)]x v. (28) 



Here and in the following the time argument t is dropped 
(because of stationarity) . The multilayer Green tensor 
constructed in terms of generalized Fresnel coefficients is 
given in Ref. 15j. The partial translational invariance 
of the problem (the layers are assumed to have infinite 
lateral extension) naturally leads to a decomposition of 
the Green tensor into an angular spectrum of s- and p- 
polarized plane waves, whose wave vectors have real com- 
ponents parallel to the multilayer surfaces. 

Let the jth layer be free space. The Casimir force 
(per unit area) between the two stacks separated by that 
layer is then determined by the zz-component of the 
stress tensor obtained from T(r, r') in the coincidence 
limit r' — > r, 



T,,,j{r,r) = lim T,^j(r,r') 



(29) 



(the index j of a quantity indicates that the quantity 
refers to the jth layer), where - as already mentioned 
- the (divergent) bulk contribution to the Green ten- 
sor must be dropped. Straightforward calculation yields 



'2h 



s,p, Af = ±, N = ±) 



dio coth 



2fcBT 



(30) 



Here, q is the transverse component of the wave vec- 
tor kj = (q, /3j), whose z-component ('propagation con- 
stant') /3j is given by 



(31) 



and 



5™^(z,z',g,w) - exp(Mi/3jz)exp(iVi/3jz') (32) 

is related to the scattering part of the Green tensor as 
[r=(p,^)] 

Gf^*(r,r',c^) = y d2ge'i(''"''')Gf^'(z,z',q,w), (33) 
Gf-*(z,z',q,..) 
= E E (-,^',9,^)e^(q) ® e^^-(-q), (34) 

a M,N 



with polarization unit vectors 



e±-(q) = ^xe„ e±-(q) = (^ge, T ft^ j (35) 

[e^ , unit vector in z-direction] . Making use of the explicit 
form of the Green tensor as given in (15| . we find after 
some algebra 



/3| E E^^5, 

i(5 



MN 
<^3 



where 



Daj = Da]{q,L0) = 1 - 



j + 'i- 



(37) 



can be thought of as accounting for multiple reflections, 
with the generalized Fresnel coefficients rjj_ = rjj_ (g, u) 
being the reflection coefhcients for cr— polarized waves at 
the top (-I-) and bottom (— ) of the jth layer. They can be 
calculated recursively (for useful recurrence relations, see 
Appendix and in this way expressed in terms of the 
thicknesses and permittivities of the layers that are actu- 
ally under consideration. Inserting Eq. (|36|l into Eq. I|3UI) 
eventually yields 



" 2^2 



duj coth 



huj 
2fcBT 



X ^ 



/•OO 



(38) 



(3t, real part). Since T^^j does not depend on the space 
point in the jth layer, the argument r has been dropped, 
and Eq. gives the Casimir force (per unit area) that 
acts on arbitrary multilayered walls. 

It should be pointed out that in Eq. nothing is said 
about the details of stratification of the walls. In fact, 
any stratified system (whose material properties change 
only along one direction - here the z-direction) admits of 
separation into an angular spectrum of s- and p-polarized 
fields which do never mix (iq . This implies that Eq. H38|l , 
where the Casimir force is expressed in terms of the re- 
flection coefficients of the walls, is very general and also 
applies to walls whose permittivity is a (on a macroscopic 
scale) continuously varying function of z. For such walls, 
however, the reflection coefficients cannot be calculated 
from simple recurrence relations. In any case, they might 
be determined experimentally. 

When in the zero-temperature limit the overall system 
is in the ground state |0), so that i{r,uj) |0) = Vr,w, 
then the thermal weighting factor coth[fiii;/(2fcBT)] does 
not appear in Eqs. H27|) and (|28|l and in the equations that 
follow from them, and thus Eq. H38|) changes to Tomas' 
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formula 



h 

'2^ 



doj X 



X 



(39) 



2. Imaginary frequencies 

Exploiting the analytical properties of the w-integrand 
in Eq. H38|l in the upper complex frequency half-plane, 
we can equivalently rewrite as 



d£. 



cothi ^ I X 
2fcBTy 



•^0 



(40) 



i.e. the frequency integration is performed on a straight 
line parallel (and infinitesimally close) to the imaginary 
frequency axis (for details, see Appendix [2)). Note that 
(the small) r] in Eq. (|40() only indicates that the (first- 
order) poles of the hyperbolic cotangent at the imaginary 
frequencies 



= i(,m ~ 2im'Kk^T /h 



(41) 



(to, integer) have to be kept to the left of the integration 
contour. 

In order to further evaluate Tzzj" as given by Eq. (|40|l . 
we first note that in the limit 77 0-1- the hyperbolic 
cotangent becomes purely imaginary. Since the permit- 
tivity is purely real and positive on the (positive!) imagi- 
nary frequency axis (see, e.g., the propagation con- 
stant becomes purely imaginary, 



(42) 



and thus the generalized reflection coefficients become 
real. Hence, the intervals between the poles (|41|l do 
not contribute to the imaginary part in Eq. (|40|l . since 
the term within the square bracket becomes purely real. 
Clearly, the same conclusions can also be drawn directly 
from Eqs. H27() and (|28|l . because the imaginary parts 
of both the permittivity and the Green tensor vanish 
at imaginary frequencies. In this way, the ^-integral in 
Eq. H40I) can be given by a residue series according to 
(r?^0+) 



with the poles from (|41() . Taking into account that 



lim {lj 



) coth 
d 



tanh 



2fcBT 



2kBT/h, (44) 



[the rest of the integrand is holomorphic; see Appendix 
we eventually derive 



7 rp ^ 

Jo 



(45) 



which may be regarded as a generalization of the popular 
Lifshitz formula (5|. Note that the zero- frequency term 
in the Lifshitz formula has been a subject of controversial 
debate, because the reflection coefficients can be discon- 
tinuous at the point {q = 0,u! = 0), so that it matters from 
which direction this point is approached |20| . From the 
above given derivation of Eq. (|45|l it follows that first the 
g-integral for a small but non-zero value of ^0 should be 
calculated and then one can let ^0 0+- It should be 
mentioned that the feasibility of flipping the contour of 
the frequency integration results from the properties of 
the Green tensor in the position space, not in the (q, z)- 
space. 

To obtain T^z.j in the zero-temperature limit, we may 
simply set (thanks to 77) T = in Eq. . The thermal 
weighting factor thus reduces to unity and we have 




dqq X 



'1+ 



>, (46) 



u—i^ 



which we may rewrite, on using Eqs. H37(l and H42|l . as 



^0 

X 



i- i+ 



(47) 



(jj—i^ 



Of course, Eq. H47() can be obtained by replacing 
Em=o • • • ^ith fi (27rfcBT)-i • • • in Eq. ggi, since 

the distance A^^ = 2iTkBT/h between neighboring poles 
[see Eq. (|41|) ] becomes small in the zero-temperature 
limit, or, alternatively, by flipping the frequency integra- 
tion contour into the imaginary axis directly in Eq. H39I) 

M. 



= 27r3 



h J2 (1- 5MRes/(o 



m— 



3. ID systems 

(43) In order to compare the above given macroscopic ap- 

proach to the Casimir force with the more microscopic 
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approach developed in j8| for one-dimensional systems, 
we have to perform our calculations for one field compo- 
nent only, e.g.. 



E(r,w) =^-i/2^(z,c^)e. 



(48) 



(A, normalization area; e^;, unit vector in x-direction). 
According to the quantization scheme outlined in Sec- 
tions III Al the (effectively) scalar electric field strength 
E_{z,Lu) can be expressed in terms of scalar bosonic ba- 
sic fields and a scalar Green function. Following the line 
in Sections III Al and III Bl and restricting, for simplicity, 
our attention to the zero-temperature limit, we derive 
(Appendix Id]| 



T 



ttA 



(49) 



with A{-^oo) being the normalization area perpendicular 
to the z-direction. Note that instead of Eq. Ij^lfl now 



(50) 



is valid, and the polarization index a can be omitted, 
since Eq. (|48|l implies normal incidence and fixed polar- 
ization. Equation (|49|l can formally be obtained from 
Eq. H39|) by making the replacement 



is 



(51) 



and rewrite Eq. (|47|l as 

F = T . - j!fL [ ,7k K^p-^'^j'ij ( 



-it/2 



a ^ ~ exp(-2Kjdj) 



(i0 sin (j) X 

(54) 



Here and in the following we do not explicitly indicate 
that rj_|_ — (w = icKj cos (j),q = Kj sin 4>) . Now, making 
the Kj-integral dimensionless by letting 



du 



-2d, e'^''''^^ dn. 



(55) 



at once extracts the (generic) asymptotic dependence on 
dj and also yields a finite domain of integration (recom- 
mended for numerical computations) 



F = -Foid,) X 
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where 



it/2 



^ 1 - r'^_rf,u 



hcTT^ 1 

240 S 



(56) 



(57) 



is the well-known formula for the force (per unit area) 
between two perfectly reflecting walls, which was first 
derived by Casimir |3|. 



and keeping only the normally incident waves (q = 0) of 
fixed (s or p) polarization a. 

In the case of two identical plates (separated by vac- 
uum, i.e. Sj — 1), the refiection coefficients Vj^ and rj_ 
can be identified with the single-plate reflection coeffi- 
cient Tj, i.e. rj± =rj. From Eq. H49(l it then follows that 
the Casimir force (per unit area) resulting from the two 
possible polarizations, Fid = '2Tzz,j, can be given by 



Fu 



2h 

TTCA 



duj 



j,2 ^2iujdj / c 



'• exp {2iujdj/c) 



(52) 



which corresponds to the result obtained in SJ, if the 
exact permittivity of the plate material is identified with 
the model permittivity in @ . Note that FioA is the total 
force acting on the normalization area A (— *■ oo) and not 
the force per unit area (as erroneously stated in 0). 



III. ASYMPTOTIC DISTANCE LAWS 

Let us consider the zero-temperature limit in more de- 
tail. In order to evaluate Eq. H47|l for the case when ej = 1 
is valid, it is appropriate to change to polar coordinates. 



Kj cos ( 



(53) 



A. Standard long-distance law 

Appendix IbI shows that the reflection coefficients rj_|_ 
as functions of q and w, rj_(_ =rj_|_(g, w), do not depend 
on the distance dj between the two reflecting walls. Any 
dependence on dj of the integral expression in Eq. H5fci|) 
is therefore a matter of how the reflection coefficients 
scale with the variable u. If dj is large, then only the 
values of rj_|_ for sufficiently large wavelengths (i.e. small 
values of both q and can effectively contribute to the 
integral expression in Eq. (|56l) . whereas the other ones 
are exponentially suppressed. 

Introducing the 'static' values of the (^-integrals of 

(r-r-_)™= lim / sin0 (rj+rj^)" , (58) 

we may evaluate Eq. H54() [or Eq. H56|l ] in the large- 
distance limit to obtain 



F 



2C(4) 



5]Li4(rJ+rJ_) (dj^oo). 



(59) 



where we have expanded 1/ [1 - 



of rj^rj. 



Tj^r^.e 



m powers 



and replaced the resulting 0-integrals 
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of powers of rj_^_rj_ by their 'static' values according 

to Eq. HHIl- In Eq. L^l^) = Em=i ^"M" is the 

polylogarithm function (the series converges for \z\ < 1, 
s > 1), and (^(s) = X)m=i is the Riemann zeta func- 
tion [C(4) = 7r''/90]. Clearly, Eq. ^ gives the correct 
asymptotics only if the 'static' polylogarithm does not 
vanish for the two polarizations (see Section lllllj(l . Us- 
ing the relation 

Li.(|x| < 1) < C(s) =Li,(x^ 1-) (s>l), (60) 

from Eq. H59() we see that the asymptotic value of the 
force is bounded by Fg {dj ) . Note that assuming constant 
reflection coefficients would formally produce the well- 
known dj'^ distance law for arbitrary dj. However, this 
unphysical assumption clearly contradicts the validity of 
Eq. jll- 

Let us briefly discuss the validity of Eq. (|59|l . In order 



We substitute this expression into Eq. (|54|l and find that 
the leading term of the force in the large-distance limit 
now reads as 



15hc ^ 



where 



= lim 

^ K,— 0+ ^ 



tt/2 



(65) 



(66) 



Thus, the Casimir force asymptotically behaves like 
^ dj^ . Obviously, other non-standard large-distance 
laws can also be observed. In particular, when the rela- 



to replace (r-J+rJ__ )" by (rj+rj_ )- according to Eq. ® , tion ^ is valid for either rf_ or rf, , then the Casimir 



the reflection coefficients rj_j_ as functions of Kj must be 
slowly varying on a Kj-scale of the order of magnitude 



of d: 



From the structure of the reflection coefficients 



(Appendix 0) it is seen that they depend on K,j via the 
(dependence on frequency of the) permittivities £i of the 
layers and the exponentials exp(— 2Ki(ii). Hence, two con- 
ditions must be satisfied. If is the characteristic fre- 
quency scale of variation of ei on the imaginary frequency 
axis, then Kj « d~^ must be small compared with ^;/c. 
Thus, one condition can be given by 



c 



(61) 



(ly^j), i.e. the characteristic wavelength dj of the 'cav- 
ity' formed by the two multilayered walls must be much 
larger than the characteristic wavelengths of all the wall 
permittivities .21.] . The other condition comes from the 
requirement that Kidi <C 1 on the relevant K^-scale men- 
tioned above. Recalling Eqs. (|42|l and (|53() and the con- 
dition (ED), we thus arrive at the condition that 



dj > V^ii^^i) di 



(62) 



{I y^Q, j,n). Note that in the case of semi-infinite walls 
considered by Lifshitz only the condition (|()T|l is 
needed. 



force asymptotically behaves like ^ d~ . 

To be more specific, let us consider the reflection co- 
efficients r|_|_ in more detail. By assuming finite values 
of ej±i(0) and dj±i, from Eq. (together with the 
properties of the permittivity) it follows that {sj ~ 1) 



< 



Kj — Kj±i 



< 0. 



(67) 



Thus, using Eq. ljB6|l (and the corresponding equation 
for r*_) and writing 



Kl + Kl±l 



'{l±l)± 



'l± 



Kl±i 2K,±id,±i s 



(68) 



we can establish that for I ~ j the inequalities 



'3± 



> 



K,±l 



'0±1)±| 



(69) 



and 



B. Non-standard long-distance laws 

Let us consider, e.g., the case when for small values of 



Kj the relation 



(63) 



is valid, with the -RJ^ being bounded functions of 0, so 
that we may write 



1 — exp{—2Kjdj) 



= K'^jR'^_R'^+ + 0(4) (64) 



r3±\<l 



Hj '*j±i 



e ^(j±i)± 



Kj + Kj±i 

are valid. The inequality H69(l shows that if 



lim \r^,+i)+ \ ^ lim 



Kj -\- Kj^\ 



(70) 



(71) 



holds, then (for the chosen values of (j)) approaches 
non-zero values in the limit Kj — >0-|-, making Eq. Ht)3|) 
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impossible. Hence, the law [Eq. can be ex- 

pected to hold for large distances. By contrast, if there 
are single-slab (dielectric) walls, then Eq. yields for 
l^j±l 



0±i)± 



Kj±l + Kj 



(72) 



(J±2)± 



:0), and the inequality H70|l thus implies that 
■|_(_ vanishes (uniformly with respect to (jj) as 0{Kjdj±i) 
in the limit Kj —>■ 0+, and a behavior as in Eq. is 
observed. Note that from Eq. H42|l for kj and the cor- 
responding equation for k^+i it follows that the relation 



Hj+i = Kj ^1+ [ej+i{i£,) - 1] cos2 (73) 

is valid, which reveals that kj+i vanishes with vanish- 
ing Kj . Recall that according to Eq. Kj 0+ en- 
tails ^^0-|-. In a similar way it can be shown that 
(for single-slab walls) the reflection coefficients for p- 
polarization, r^±, also vanish uniformly as 0{Kjdj±i) in 
the limit Kj— >0-|-. Consequently, when the two walls 
are single-slab dielectrics, than the rj_|_ behave according 

to Eq. (jnSl, and hence the d^^^-law is observed for large 
distances, with the functions Rj^{4>) in Eq. H65|l being 
proportional to the respective slab thickness. Clearly, if 
only one of the two walls consists of a single slab, then 
the dj^-\aw may be observed. It should be pointed out 
that the large-distance asymptotic regime again requires 
the conditions and Ht)2|l to be satisfied. 

Let us remark that it is conceivable that other special 
choices of the walls may produce other than ~ dj^, 

dj^, and ^ dj^ asymptotic distance dependences of the 
Casimir force. It may also happen that in the asymptotic 
expansion of the Casimir force terms d"" with different 
values of n must be taken into account for not extremely 
large distances, if the weights of the terms substantially 
differ from each other. 

We finally note that the one-dimensional counterpart 
of the dj'^-law is a d~^-\sm. It changes to a dJ'^-lsLW when 
rj± ~ £,Rj±/c holds in the limit ^— >0-|-, and it changes 
to a d~^-\a,w when only one of the reflection coefficients 
shows this behavior. In particular, two single-slab walls 
give rise to a d~^-law in place to the standard d~^-\a,w. 



may be approximately replaced by single-interface reflec- 
tion coefficients according to 



'3± 



Kj/Kj±l ~ l/Sj 



±1 



Kj/Kj±i + l/ej±i 



(75) 



Thus, we effectively deal with two semi-infinite walls of 
permittivities Sj±i, so that Lifshitz' approximation for 
short distances can be used |^: 



Kj±l 



and 



'j± 



0, 



'3± 



£,-±1 + 1 



(76) 



(77) 



Note that Eqs. H74|) and H76|) are consistent with each 
other, and Eq. (|76() is valid if 



ej±i(cj = iKjc) — 1 <C 1, 



I.e. 



(78) 



(79) 



where the plasma frequencies are defined by 

flf = lim;j^oo(aj^ [1 — ei{uj)]) 18]. In this approximation, 
Eq. H47|) takes the well-known form of {v = 2qdj) 



167r2df JO 



X / dvv 
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d^ X 

(£j + l + l) (£j-l + l) ^ 

{£]+! - 1) - 1) 



(80) 



Let us consider, for simplicity, single-resonance media of 
Drude-Lorentz type, such that 



(81) 



For small 70, Eq. H80|) can then be further evaluated to 
obtain (Appendix El 



F : 



172/2 Lis 



(82) 



with 



C. Short-distance law 

For short distances, we have to compare the Kj-scale of 
variation of the reflection coefficients with the now 

large scale of variation dj^ of e"^"^''^. In particular, 
assuming dj <^ dj±i, so that we may let 

from Eqs. (jB6p and ljB7|l (and the corresponding equa- 
tions for rf_) we see that the reflection coefficients rj_|_ 



~. 1 ^ r(4m-l) £ 

2 ^ [r(2m)]2 m3 



Z-^- (83) 



IV. NUMERICAL RESULTS 

In order to illustrate the dependence of the Casimir 
force on the various parameters, we have evaluated 
Eq. (|56|l numerically for single-slab walls fFigs. El- I15|) 
and a periodic multilayer wall structure (Figs. 1161 and 
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10 

10~^ 10° 10^ 

FIG. 2: Dependence of the relative Casimir force F/ Fq be- 
tween two identical single-slab walls on the thickness d± of 
the walls and the wall separation d (material parameters: 
LOO = 1.0 X 10^ s-^ = Lui^o, ^ = 1.6176 X 10^'' s"^ ^ IO'^wlo, 
70 = 9.7 X lO^'' s~^ ^ l(f ojLo). 

I17|l . The results extend, in a sense, the ID results given 
in 1^ to three dimensions. In Figs.[21-El the wall mate- 
rial is characterized by a single-resonance Drude-Lorentz 
permittivity of the type given in Eq. H81(l . We have per- 
formed the calculations for transverse resonance frequen- 
cies (jjQ in two qualitatively different frequency domains, 
namely 

Wo =i WLO = 1-0 X lO^s"^ (84) 

and 

tjo =it^Hi = 2.0 X lO^^s-^ (85) 

While the higher frequency whi corresponds to dielec- 
tric material (such as Si), the lower frequency wlo may 
be regarded as being typical of metal-like material (such 
as Mg). To allow a comparison with the ID results in 
Rcf. |9j, we have performed the numerical calculations 
for the parameter values used therein. Note that in the 
contour plots there are always 15 lines of equal relative 
force F/Fq, equidistantly between the highest and the 
lowest occurring value. 



10-2 

10-^ 10° 10^ 

d± (Mm) 

FIG. 3: The same as in Fig.H but with ujo = 2.0 x lO^'^ s"^ = 
c^Hi, n = 6.536x10^^ s"^ « 3wHi,and7o = 9.859x10^^ s"^ « 
0.01 a;Hi. 



Let us first consider the relative force F/Fq between 
two single-slab walls. The influence of the wall thickness 
d± =dj±i on the distance law for the two resonance fre- 
quencies Wo —i-^ho and luq =ljhi, respectively, is shown in 
Figs. Eland 13 (d = rfj), and the dependence on ujq of F/Fq 
for chosen d± is shown in Figs. 0] and [3 From Figs. |2| 
and O it is seen that (for sufficiently high resonance fre- 
quencies) the relative force increases with the distance 
between the plates, attains a maximum, and eventually 
decreases with further increasing plate separation. With 
increasing thickness of the plates the maximum becomes 
broader and a plateau-like behavior is observed. Fig- 
ures El and 01 reveal that for low resonance frequencies 
the plateau can become very broad, so that observation 
of decreasing values of F/ Fq would require very large dis- 
tances, at which the force effectively vanishes. Compar- 
ing Figs. El and O we see that the response of the force to 
a change of the plate thickness is much more sensitive for 
high resonance frequencies than for low ones. Needless 
to say that for sufficiently thick plates the force becomes 
independent of the plate thickness. From Fig. it is 
seen that the long-distance asymptotic behavior F ^ 
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10" 10^ 10-^ io» 

FIG. 4: Dependence of the relative Casimir force F/ Fq be- 
tween two identical single-slab walls on the resonance fre- 
quency ojo and the wall separation d (wall thickness d± — 



0.5 fim; material parameters: 
lO^cJLO, 70 = 9.7 X 10^'' s"^ f» 



10^ 



1.6176 X 10'" s" 



[Section IlII B| is observed when - in agreement with the 
condition H()2|) - the distance between the plates substan- 
tially exceeds the plate thickness. A comparison of the 
3D results in Figs.[21and|21with the ID results in Fig. 2 in 
Ref. jOl] shows a quantitatively rather than qualitatively 
different behavior of the relative force F/Fq in the two 
theories. Clearly, the force itself behaves quite different 
in the two theories. 

Figures 0] and |31 clearly show that increasing the reso- 
nance frequency loq generally lowers the force. In partic- 
ular, it is seen that the position of the maximum of the 
(relative) force is shifted to smaller values of the wall sep- 
aration when the resonance frequency increases. At the 
same time, the maximum value decreases and the long- 
distance asymptotic behavior sets in at smaller distances. 
The dependence of the force on the resonance frequency 
is in agreement with the Drude-Lorentz permittivity (|81|l . 
For chosen plasma frequency Q and absorption parame- 
ter 7o, the maximum absolute value of the permittivity 




FIG. 5: 



0.5 1.0 l.i 



The same as in Fig. |1] but with d± 
« Sum, and 70 = 9.859x10^^ I 



2.5 ^m, Q — 
^ ^ 0.01 tJm. 



decreases with increasing value of loq, thus reducing the 
reflection coefRcients of the walls. From Figs. and [3 
it is seen that increasing the value of the resonance fre- 
quency has a similar effect as decreasing the value of the 
wall thickness. Both a very small plate thickness and a 
low permittivity can lead to poor plate reflectivity. 

The influence of the plasma frequency on the distance 
law for chosen plate thickness and the two resonance fre- 
quencies ujo =wlo and loq =tom is shown in Figs. |H1 and El 
respectively. Figures [TUIandlTTIillustrate the dependence 
of the force on the plasma frequency and the thickness of 
the plates for chosen distance between them. Since the 
plasma frequency can be regarded as being a measure 
of the strength of the medium resonance, higher values 
of the plasma frequency imply higher values of the plate 
reflectivity and thus higher values of the force, as can 
be clearly seen from the figures. Moreover, the width of 
the band-gap featured by the permittivity H81I) increases 
with the plasma frequency, which explains that the in- 
terval of large (relative) force also grows with increasing 
plasma frequency (Fig.O. In particular, Figs. 1101 andllll 
show that a high plasma frequency can compensate for a 
small plate thickness. Note that ^ oo corresponds to 
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FIG. 6: Dependence of the relative Casimir force F/ Fo be- 
tween two identical single-slab walls on the resonance fre- 
quency ujo and the wall thickness d± (wall separation d — 
1 /im; material parameters: Q, = 1.6176 x 10^^ s~^ ~ 10* <j-'lo, 
70 = 9.7 X 10" s"^ « 10*^ a;Lo). 

perfectly reflecting plates, whose thickness can then be 
arbitrarily small. 

The effect of material absorption is illustrated in 
Figs.1121-1151 The dependence of the distance law on the 
absorption parameter 70 for chosen plate thickness and 
the two resonance frequencies = wlo and loq = whi , re- 
spectively, is shown in Figs. 1121 and ll3l and Figs. 1141 and 
[T51 present the dependence of the force on the absorp- 
tion parameter and the thickness of the plates for chosen 
distance between them. As expected, the force is seen 
to decrease with increasing absorption parameter. The 
effect is similar to that observed when the resonance fre- 
quency is increased. In both cases the maximum absolute 
value of the permittivity decreases (for chosen plasma 
frequency), so that the reflection coefficients of the walls 
diminish. It is worth noting that the force responds more 
sensitively to a change of jq/luq for high resonance fre- 
quencies than for low ones. In particular, from Fig.ll3lit 
is seen that in the first case the force is practically not 
influenced by material absorption as long as 70 /luq «C 1 is 
valid, and it effectively reduces to zero when 70/wo sub- 
stantially exceeds unity. Figure El also shows that the 




5X10-'' 5x10-2 5.0 

Wo/WHI 



FIG. 7: The same as in Fig. [S] but with d = 1 /im, Q — 
6.536x10^^ s-i «3d^Hi,and7o = 9.859x10^^ s"^ wO.OIojhi. 
The lowest contours are not trustworthy, because of numerical 
errors. 



maximum of the relative force F/Fq is shifted to larger 
values of the distance between the plates when the value 
of 7o/a;o increases. For low resonance frequencies the 
ratio 70/(^0 must increase to rather extreme values be- 
fore the force substantially diminishes, as it is seen from 
Fig. El Note that for a metal {ujo ^0), 70 becomes in- 
versely proportional to the conductivity. Small loq and 
small 7o, i.e. high conductivity, lead to such a broad 
plateau of nearly constant (maximum) value of F/Fq w 1 
that effectively Casimir's formula H57|l applies. 

Let us finally consider the force between two identical 
multilayered walls composed of identical bilayers, where 
each bilayer is made of a metal-like and a dielectric-like 
material (Figs. [TBI and I17() . For comparison with the ID 
results in Ref . 9] , we have again performed the numerical 
calculations for the parameter values used therein. From 
Fig. El it is seen that for sufficiently small distances be- 
tween the walls the force does not depend on the number 
of bilayers the walls are composed of. In this case, only 
the inner layers essentially determine the force, which is 
in full agreement with Eq. H75|l . 
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10^ 10^ 10^ 10* 10^ 



FIG. 8: Dependence of the relative Casimir force F/Fq be- 
tween two identical single-slab walls on the plasma frequency 
Q. and the wall separation d (wall thickness d± = 0.5 ^m; 
material parameters: loq = 1.0 x 10^ s~^ = cjlo, 7o = 
9.7 X lO^^'s-i lO'^iJLo). 



For larger distances between the walls the relative force 
increases with the number of bilayers. The changes 
in the curvature of the curves in the figures indicate 
that (for chosen number of bilayers) the distance law 
can drastically change several times before the relative 
forces becomes constant. These changes, which are less 
pronounced when the inner layers are metal-like ones 
[Figs. [TBT al and ITTT a)] . may be regarded as being typ- 
ical of a 3D theory (cf. Ref. 0)- 



V. SUMMARY 

We have studied the Casimir force between dispers- 
ing and absorbing multilayered dielectric plates. On the 
basis of the quantization scheme for the electromagnetic 
field in causal media as given in Ref. we have ex- 
tended the recently derived zero-temperature result 
to finite temperatures. The derived formula generalizes 
Lifshitz' formula Q to arbitrary multilayered walls, and 




FIG. 9: The same as in Fig. |H1 but with d± — 2.5 /im, uo ~ 
2.0 X 10^^ s~^ = LJHi, and 70 = 9.859 x 10^^ s"^ f» 0.01 a;Hi. 



application of the ID version of the theory to single-slab 
walls yields the results in Ref. 8] . 

Restricting our attention to the zero-temperature li- 
mit, we have studied the problem of asymptotic distance 
laws of the Casimir force. We have shown that Lifshitz' 
approximation for short distances also applies to multi- 
layered walls and, depending on the wall structure, the 
distance law can drastically change with increasing wall 
separation. In particular, for two single-slab walls the 
Casimir force tends to behave like d^^ instead of as 
the wall separation d goes to infinity. 

Assuming permittivities of Drude-Lorentz type, we 
have finally presented a number of numerical results in or- 
der to illustrate the dependence of the zero-temperature 
Casimir force on various system parameters and to com- 
pare with the ID results recently reported in Ref. |^. 
While in the case of single-slab walls the 3D and the 
ID theory yield a qualitatively similar dependence of the 
relative Casimir force on the wall separation, significant 
differences (not only for the absolute force but even for 
the relative force) may be observed for more complicated, 
multilayered walls. 
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5x10° 
5x10-1 fll 




FIG. 10: Dependence of the relative Casimir force F/Fo 
between two identical single-slab walls on the plasma fre- 
quency n and the wall thickness d± (wall separation d — 
1 jj,m; material parameters: loq = 1.0 x 10® s~^ — ojlo, 
70 = 9.7 X 10" s"^ ^ 10^ a;Lo). 
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APPENDIX A: DERIVATION OF EQS. 



It is convenient to discretize the fundamental field vari- 
ables f (r, uj) so as to work with H — ^ hio^ ^^fR./j 
instead of Eqs. (|19|l [(R, w^), grid points], and hence with 
a partition function of the form Z = Y\^^ ^ ^. From 
the continuum limit of such a calculation, it then easily 
follows that 



(ft(r,o.)®f(r',c.')) 



coth 



2kBT 



- 1 




5x10° 5X10^ 5x10^ 

O/WHI 

FIG. 11: The same as in Fig. 1101 but with d — 1 /im, ujo 
2.0 X 10^^ s"^ = tJHi, and 70 = 9.859 x 10^^ s"^ ^ 0.01 wm. 



(f(r,w)®f+(r',cj')> 



coth 



2kBT 



(A2) 



(Al) 



as well as 

{i{r,uj)(g)i{r',Lj')) ^{i^r,ou)(g)i\r\Lu')) ^0. (A3) 

Substituting Eq. and the corresponding equation for 
the displacement field into Eq. H23I) . we may write, on 
recalling Eqs. (221), (CH), dl) and |J3)l 

ri(r,r',t) 

= / duj dw' e-*('^-'^')*(D(r,w)(g)Et(r,t^)) 
Jo Jo 

-f e'('^~'^')*(Dt(r,u;)(g)E(r,u;))j . (A4) 

Making here explicitly use of Eqs. 1)13(1 and p4|l and the 
reciprocity property we obtain 

ri(r,r',i) 

/■OO /"OO 

— fJ,o doj I div' uj'^V X V X K{r,r' ,oj,oj' ,t), 
Jo Jo 

(A5) 
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To/l^LO 

FIG. 12: Dependence of the relative Casimir force F/ Fo be- 
tween two identical single-slab walls on the absorption pa- 
rameter 70 and the wall separation d (wall thickness d± = 
0.5 /xm; material parameters: = 1.0 x 10® s~^ — uj-lo, 
n = 1.6176 X 10^*^ s"^ « 10* WLo). 



where the tensor-valued function If (r, r', u), lo' , t) is given 
by 

i<:(r,r',w,w',t) 

= e~'i^-^')\{ G ★ VKf ) (r, u) (g> (Wf^ ★ G*) (r', lo')) 

+ e*(---')*( ( G* * W^f (r, c^) ® {Wt ★ G) (r', cj'))- 

(A6) 



Here, the abbreviating notations 



(A7) 



and 



(G*f)(r,w) = y d=^r' G(r,r',a;)f(r',w) (A8) 

have been used. Applying the integral relation (|17|l and 
using Eqs. (jAl|) and HA2p . from Eqs. (jA5|) and (|A6|) we 





6X10-** 5X10—'' 5X10° 
7o/i^HI 

FIG. 13: The same as in Fig. [HI but with d± = 2.5 fim, 
ujo = 2.0 X 10^^ s"^ = um, and = 6.536 x 10^^ s~^ ^ Sum- 



derive 
ri(r,r') 

= ^l^du; coth(^^) V X V X 3[G(r, r', 

(A9) 

Recalling Eq. p5|l . we see that Eq. (|A9|I just leads to 
Eq. The magnetic part T2(r,r',i) [Eq. ^] can 

be calculated analogously. In place of Eq. (|A4|) we now 

have 



= dw dcu' e-^('"-"')*(B(r,w)®Ht(r,c^)) 
Jo Jo '- 

+ e'("-"')*(Bt(r, u) ® H(r, cj))] , (AlO) 
from which by means of Eqs. H12|l and (|14|l follows that 
T2(r,r',i) 

= — /zo / dui duj' ujuj'V X K{r,r' ,uj,iv' ,t) X V' . 
Jo Jo 

(All) 
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7o/l^LO 

FIG. 14: Dependence of the relative Casimir force F/Fo 
between two identical single-slab walls on the absorption 
parameter 70 and the wall thickness d± (wall separation 
d — 10 /im; material parameters: uq — 1.0 x 10^ s~^ — uj-lo, 
Q. = 1.6176 X 10^^ s"^ « 10** LULo)- 

Comparing Eq. (|A11|) with Eqs. (jA5|) and HA9|I. we see 

that 

r2(r,r') 

t pOC / Pi \ < 

^--J^ dc.coth(^^j Vx3[G'(r,r',c.)]xV' 

(A12) 

which is just Eq. (^5)) . 

APPENDIX B: REFLECTION COEFFICIENTS 

Let Ef{z) be the classical (complex) electric field that 
is observed in the Zth layer when in a layer with I' <l an 
electromagnetic wave (of chosen frequency tu, transverse 
wave- vector component q and polarization a) propagates 
in z-direction (i.e. upwards in Fig.^). It can be written 
as 

EUz) = El,t^,^i [e+ e^A^ + e^'^^Tr+e" e^^^] , (Bl) 




To/l^HI 

FIG. 15: The same as in Fig. ^] but with d = 10 fim, luq = 
2.0 X 10^^ s"^ = WHi, and Q = 6.536 x 10^^ s"^ ^ 3ujm- 

where is the amplitude of the upwards propagating 
wave just before the upper boundary of the I'th layer, 
t^i^i is the transmission coefficient from the I'th to the 
Zth layer {I' <l), and rf^ is the reflection coefficient at 
the upper boundary of the Ith layer (e^j according to 
Eg. a = s,p). Note that we adopt the convention 

[l5l | that z = denotes the lower boundary in all layers I 
with I > 0, and z = di denotes the upper boundary in all 
layers I with / < n. The transmission and reflection coef- 
ficients are determined by the requirement that etEf{z) 

and (et x V)i?^(z) are continuous at the surfaces of dis- 
continuity (ej, tangential unit vector). Applying V in 
the (q, z)-space as iq + Bzd/dz, straightforward calcula- 
tion yields the sets of equations 

i?,/,e^^''^'[l + rf+] 

= Pi+it', -e''^'-^'''^%+,^+] (B3) 
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FIG. 16: Dependence of the relative Casimir force F/Fo be- 
tween two identical multilayered walls on the wall separation 
d for various numbers n of identical bilayers (a) L-H and (b) 
H-L each wall is composed of. The data of the individual 
layers L and H are respectively d± = 0.01 fim, ljq = wlo = 



1.0 X 10^ s-\n = 1.6176 X 10'" s" 
d± = 2 fim, uo = LJm = 2.0 x lO^'^ ; 
70 = 9.859 X lO^^s^V 



70 



9.7 X 10" s"^ and 




FIG. 17: Detail of Fig.lTCl 



[fci=a;e,Hi/2/c=(g2+/32)i/2] and 



[1 



ki+i 



1 



(i + l) + J' 



(B4) 



for s and p polarization, respectively. Eliminating in 



Eqs. (|B2l) - (|B3) ip/je^^'^'Vip/i+i, we derive the recur- 



{(3i/(3i+i - 1) + + l)e2^A+'^'+'rfj^, 



and 



('+!) + 



(A//3/+1 + 1) + (/3///?i+i - l)e2^A+i'^'+ir«,^^^^ 

(B6) 

+ 1 ei + i> ^ '^ft + i ^ ei + i)*^ '(i+l) + 

^/3, + i 



£1 + 1'' ^ + 



(B7) 



which terminate at l = n—l, because of r5^_|_=0. Note 



that the 



(«-!) + 



for the last surface of discontinuity are 



the well-known single-interface coefficients. 

By considering in the same way the solution that 
emerges in the /th layer from a downwards propagating 
wave in the I'th layer {I' > /), analogous recurrences for 
the coefficients rf_ are derived. They are formally ob- 
tained from Eqs. I|B6I) and (|B7I) by making the replace- 
ments 



(B8) 



and turning the subscript + of the reflection coefficients 
into — . Because of rQ_ = 0, they terminate at Z = 1. From 
the recurrences it follows that the reflection coefficients 
are real on the (positive!) imaginary frequency axis and 
fulfill there the inequality 



-1 < rf± < 1 



(B9) 



It should be noted that it is possible to give the i^^^ 
explicitly once all the rf_j_ have been computed. 



APPENDIX C: INTEGRATION ALONG THE 
IMAGINARY FREQUENCY AXIS 

From Eq. (|15|1 the (exact) Lippmann-Schwinger-type 
integral equation 

G''^"*(r,r',cj) 

= ^ y d^s G^""^(r, s, w)fc(s, iu) G(s, r', u), (CI) 

can be derived, where G*"""^ and C"'''' are the bulk and 
scattering parts of the full Green tensor G, respectively, 
and Se{r, w) is the associated deviation of the permittiv- 
ity from the 'bulk' situation. Since any possible Green 
tensor (such as G or Q^^^^^j has an asymptotic behavior 
(c^/w^) for \uj\ — > oo in the upper complex half-plane 
Eq. lEU reveals that u;2e(r, w) G"'=''*(r, r', w)/c2 has 
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exactly the same asymptotic behavior as the permittiv- 
ity difference 5e{r,uj). We can therefore conclude that 
frequency integrals of type 



(C2) 



which appear in the zero-temperature versions of 
Eqs. (|27|l and (|28|l . are convergent along any contour C 
running from the origin (for metals the origin should be 
excluded, because of the pole at that point) to infinity in 
a chosen direction of the upper half-plane and yield al- 
ways the same value, provided that Se{r,uj) approaches 
zero at least as 0(l/|aj|-'^+'^) ((5 > 0) when \uj\ goes to in- 
finity, which is the case. 

Hence, the frequency integration (for zero tempera- 
ture) in Eqs. ^ and ^ (G ^^ g-scat^ performed 
along the positive imaginary axis (instead of the positive 
real axis) and the imaginary part can be taken after the 
integration has been performed. In this way, Eq. H39|) 
becomes equivalent to Eq. 14t)|) . To include the thermal 
weighting factor (which behaves like uj^^ for s-0), we 
first note that integrands of the type which appear in in 
Eqs. (EH) and (EHl (G^ C''*) remain perfectly regular 
at w = 0, thus 



= lim 



duj coth 



2kBT) 



(C3) 



Since coth[/ia;/(2fcBT)] is holomorphic and bounded for 
^ now may change the integration path in 

Eqs. (|?7jl and if^ according to 



I auj coth — — — 
Jo V2fcBT 



^ ' duj coth I 



2fcBT 



, (C4) 



which then leads to Eq. 



APPENDIX D: DERIVATION OF EQ. l|l9|l 

The ID electric field strength E(r, lS) as given in 
Eq. H48|) implies the equations 



B(r,c.) 



By dE_{z,uj) 



jy/A dz 



(Dl) 



D r, u;) = fj= -\' ' p2 

in place of Eqs. ifT^ and ifT^ . so that Eq. if^ reduces 
to 



r,,(r, r', t) = -i(D(r, t)E(r', t) + B(r, t)H(r', t)). 



(D3) 



We now proceed further as in the 3D case and substitute 
the ID version of Eq. into Eq. HD3|I . where the ID 
Green function solves the inhomogeneous wave equation 



— ^ + —e{z,uj) 



G{z,z\uj) = -d{z- z'). (D4) 



Restricting our attention to the zero-temperature limit, 
we derive, on applying the ID versions of Eqs. I|16|) and 

(0|D(r,t^)Et(r',t^')|0) 



A dz-^ 



{0\{G*Wf){z,Lo)iWfUG*)iz',Lo')\0) 



-:^S{u:-u;')^^-s[G{z,z',u;)] 

h up' 

—^8{uj - J) — '^ [£(z, lS)G{z, z', w)] 



(D5) 



and 

(0|B(r,a;)H^(r',w')|0) 



-{Q\{G^Wf ){z,uj) {Wp *G*){z',J)\ 0) 



A dzdz 

= _^5(^_^')^!iJ(G^£"G*)(z,z',c.) 
■kA ozoz' &■ 

^ 5(c.-c.')^3[G(z,z',c.)] 
azaz' 



ttA 



(D6) 



[for W and the ★-notation, see Eqs. HA7|) and (jA8|) ]. By 
Fourier transforming the fields back into the time domain 
and combining with Eq. ljD3p . we arrive at 



(D7) 



T,,(r,r') = --^ / da; 

271 A Jo 



'^e{z, uj)G{z, z', u) + Q^Q^, C!{z, z', uj) 



Note that Eq. I|D7|I is valid for an arbitrary e{z,uj). 

Let us now specify e{z,uj) according to Eq. H2()|l . The 
bulk Green function in the jth layer reads 

Gf"^(z,z',c.)--^e*ftl^-^'l 

= -^[e''3^("-^')6l(z-z') + e-*''^(^-^')6l(z'-z)] (D8) 

[0{z), Heaviside step function], where the propagation 
constant 



(D9) 
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coincides with the full wave number. The scattering part 
of the Green function in the jth layer 

Gf z', Lu) = + <?7e-**^ (DIO) 

is the solution of the homogeneous wave equation, where 
the matching conditions 



g-e-^f"'"' = r,+e'ft'^^- (^f + + .9+) {z = d, > z') 



(Dll) 
(D12) 



must be satisfied. Thus 



(0)± 



(D13) 
(D14) 



The ID refiection coefficients rj± can be constructed as 
outlined in Appendix Obviously, they are identical 
to the 3D ones (for s-polarization) taken at q = [23| ■ 
Combining Eqs. (|D10|I and (jDlSp . we find that, on taking 

gf'^^ from Eq. JD8|| . 



(D15) 



From Eq. ljD15(l it then follows that 



w2 , , d^Gf''\z,z',uj) 

-eiz,.)Griz,z'..)+ ^a.^; 



cos[/3j(z-z')]- (D16) 



Substituting this expression into Eq. IjDTjl [G{z, z',u}) i— > 
Qscat^^^ z', uj)] and setting z' — z, we eventually arrive at 
Eq. 



APPENDIX E: DERIVATION OF EQ. (15^ 

We first note that Eq. H8UI) can be rewritten as 



F'. 



1)(£,_1-1) 



L(ej+i + l)(e,_i + l) 



From Eq. l(HT|l it follows that 



e,±i(*e) + l l + 2[e(C + 7o)+^o1/^^' 



(El) 



(E2) 



Substituting this expression into Eq. ljEl|) . we obtain 

[y=(e + 7o/2)/r!] 



F : 



Ml 



(E3) 



where 



and 



/■oo 

= / [l + 2(y2 + a2)] 

^ TT7T 



-2m 



2 ^ - 7o/4 

02 • 



(E4) 



(E5) 



Under the condition that 70 is sufficiently small, so that 
70 ^ 2il and 70 /4 < cjq , the integral (|E4p is approxi- 
mately given by 



dy 



2^™+^ 7_oo (y2 + a2 + 1/2)2 



(E6) 



which can be easily evaluated by employing the residue 
theorem. For this purpose we write 



1 



(y2 + a2 + 1/2) 
2m 



2m 



1 



{y -y+){y -y-Y'"^ ^ 
1 



(y - y-)iy - y+) 



2m -1 



where 



2/± = i^v/a^ + 1/2. 



We thus derive 



26m-i [r(2m)]2 



(E7) 



(E8) 



(E9) 



Combining Eqs. HE3() and HE9|I eventually yields Eq. 
Because of 

/■7o/2n 



E4 / [l + 2(y^ + a^)] 

^ m-^ Jo 



<^Li3[(a2 + l/2)-V4], (ElO) 

the small relative error that results from replacing the 
lower limit of integration in Eq. (|E4(I by zero can be es- 
timated according to 



where [see Eq. 



1 Li3(x-V4) 



8TTy^ Li2(x-V64) ' 



(Ell) 



(E12) 



Note that the factor f{a'^ + 1/2) is at most of order of 
unity. 



J m=l 
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